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Abstract: We study the abelian avalanche model, an analogue of the abelian sandpile 
model with continuous heights, which allows for arbitrary small values of dissipation. 
We prove that for non-zero dissipation, the infinite volume limit of the stationary 
measures of the abelian avalanche model exists and can be obtained via a weighted 
spanning tree measure. Moreover we obtain exponential decay of spatial covariances 
of local observables in the non-zero dissipation regime. We then study the zero dissi- 
pation limit and prove that the self-organized critical model is recovered, both for the 
stationary measures and for the dynamics. 
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1 Introduction 

The model studied in this paper is a continuous height version of the abelian sandpile 
model of Bak, Tang and Wiesenfeld [3]. The reason for our interest in continuous 
heights is that this allows one to add arbitrarily small amount of bulk dissipation 
to the model. The presence of dissipation makes the model non-critical, [TB], and 
therefore easier to handle. In the physics literature, one recovers the self-organized 
critical model in the limit of zero dissipation, see e.g. [13]. However, so far no rigorous 
proof of this fact has been given. 

Our setting will be the abelian avalanche model investigated by Gabrielov [6]. 
That paper focuses mainly on a deterministic dynamics. However, as pointed out 
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there, many of the basic features, such as the abehan property and the burning 
algorithm, remain the same for a stochastic dynamics. Our models will live on Z*^, 
and for simplicity, we only consider a nearest neighbor toppling matrix. It is not 
difficult to extend our results to more general toppling matrices. 

Due to the toppling rule, the continuous model still retains a lot of the features of 
the discrete model. In fact, continuous heights are only used to define the dynamics, 
and the stationary measure has a natural description in terms of discrete height 
variables. The discrete effects have been observed numerically in [2]. Our main result 
is that the induced (non-Markovian) process on the discrete height variables converges 
to the standard sandpile process as the dissipation vanishes. 

1.1 Toppling matrices, stabilization, toppling numbers 

For x,y G Z*^, we write \x — y\ = \\x — y\\i := Ylf^i \xi — yi\ and we denote x ~ ?/ if 
X and y are neighbors, that is |x — ?/| = 1. Each site x will carry a height variable 
with value in the interval [0, 2(i + 7), where 7 > is a parameter. On toppling, a site 
will give height 1 to each of its neighbors, and lose height 2d + 7, that is, an amount 
7 of height is dissipated on each toppling. This can be summarized in the toppling 
matrix: 

{2(i + 7 ii X = y; 
-1 ifx-y; 
otherwise. 

For A C Z'^, we write A^^^ = {A^J)^^y(.A 

for the toppling matrix restricted to A. 
We define the sets of height configurations in Z"^ and in A by 

X = [0, oof and Xa = [0, 00)^. 

A site X is called 7-stabie in configuration rj, if rjx < A'^'J = 2(i + 7 (otherwise, x is 
7- unstable in configuration 77). The sets of 7-stable configurations are denoted 

fi(^) = [0, 2d + 7)^' and fii^^ = [0, 2d + 7)^. 

We will also use 

= {ri e X -.ri.^ <2d + -f, X E A}. 
A toppling of site a; in A is defined by the operator T^l, via the formula 

(^S^) =^.-Ai;\ 2/GA,r/GA'A. (1.1) 

We use the same formula to define T^lrj, when rj E X. In this case, note that T^lrj 
does not alter rj outside A. 

A toppling is called 'j- legal if the toppled site x is 7-unstable before toppling. A 
sequence (xi, . . . ,x„) is called {A, ■y)- stabilizing for rj E X, ii 



2 



i) ^lyi is a 7-legal toppling of T|y]^_^ o ■ ■ ■ o T|y]^r/, l<k< 



(7) 



(ii) the final configuration is in Q 

Note that for all t] E X and for all finite A C Z,'^, a (A, 7)-stabilizing sequence 
exists. The number of times a site topples does not depend on the (A, 7)-stabilizing 
sequence, and hence there is a well-defined stabilization map , see e.g. 

[5], or P, Appendix B] for a proof. 

The result of the stabilization is related to the original configuration by 

(4'^^). = Vy- (Ai"^Ari>)r^),, y G A. (1.2) 

where N^^i] is the vector consisting of the toppling numbers associated to the (A,7)- 
stabilization of 77, i.e., {N^^ri)^ denotes the number of times x G A topples during 
stabilization of rj in A. Stabilization in volume A can of course also be viewed as a 
map 4''^ -.XA-^n^j^K 

The addition operators are defined by 

: 4''^ - ^i"^^ = S^^\V + S.), (1.3) 

where Sx denotes the vector having entry equal to one at site x and zero elsewhere. 
The addition operators commute, that is, for all x, ?/ G A, 77 G 



a 



(7) /^(7) ^\ _ ^(7) 

Ay 



This follows from the fact that stabilization is well-defined: indeed both expressions 
in (11.41) are equal to S^\r] + 6x + Sy). 

We endow X, Aa, Q, Qj>^ with the product metric 

dist(r7i,r72) = 5^2-Nmin{|(r7i),-(r72).|,l}, (1.5) 

X 

where the sum is over Z'^ or over A. 

Given a function (f : A —* (0,oo), we define a jump Markov process on stable 
configurations. The action of the generator on Borel measurable functions / : — >■ 
M is given by 

= crm = E ^(^) [/ - /(^)] • (i-^) 

The above process is described in words as follows: at each site x G A we have a Pois- 
son process with intensity (p{x) (at different sites these processes are independent). 
At the event times of this Poisson process we apply the addition operator to the 
configuration. 
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1.2 Stationary measure, Dhar's formula 

As in the discrete case, the notion of allowed configuration can be defined. A 7- 
forbidden subconfiguration (7-FSC) is a pair (W, rjw) where C Z*^ is finite, rji^ G 
Xw, such that for all y G W, 

Remark 1. Note that since the off-diagonal elements of the toppling matrix A^"'^ 
are not depending on 7, the right hand side of inequality (11. 7p is independent of 'y. 
Therefore we have the same forbidden subconfigurations for any value of So from 
now on we use the words FSC rather than 'y-FSC, and allowed rather than 'j-allowed. 

A configuration 77 G fij^'' (respectively 77 G Q'^'^^) is called 

allowed if there does not exist finite W A (respectively W C Z*^) such that the 
pair consisting of W and the restriction of rj to W is an FSC. 

Let 

7^^^) = {t] e : 7] is allowed}, 
7^(^) = {r]e fi(^) : 7] is allowed}, 

= {ve n^^) : T]v G 7^S7^ for all finite V CZ'^}. 

The results of [HI Sections 3,4] imply the following properties of allowed configu- 
rations. 

Proposition 1. (i) a^^, maps T^^"* one-to-one and onto itself. 
(11) Volume (7^^^) = det(Ai^^). 

(Hi) Lebesgue measure on TZ^^^ is invariant under a^l\, x G A. 
Hence, the probability measure on n^2^ defined by 

mJ'^(A) := 



Volume (7^S■^^' 



-A J 

is stationary for the Markov process defined in Section II. 1[ 

For rj G ^^/^\ we denote by n^^\x,y,ri) the number of topplings at y needed to 
stabilize rj -\- S^- Then due to (11.21) we have the relation 

(a^^)^ = Vz + {S.)z - 5^ AWni'')(x,y,r?). (1.8) 

Averaging with respect to mj^^ and using stationarity under the action of the addition 
operators (Proposition [TJ^iii)) gives the "Dhar's formula" [H]: 

E^(.)(ni'')(x,i/,r^)) = (Aj));^V|Vn =: G^)(x,y) (1.9) 
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If 7 > or (i > 3, the inverse G^'^\x, y) := also exists and is equal to the limit 

lim^l^d y). Note that G^'^\x,y) equals the Green's function of a continuous 

time random walk that crosses an edge at rate 1, and is killed at rate 7. Similarly, 
G^^\x,y) equals the Green's function of a continuous time random walk that crosses 
an edge at rate 1, is killed upon exiting A and is killed (inside A) at rate 7. 

Markov's inequality and (11.91) imply 



m 



A 



^'^(n^;^\x,y,v)>l)<G^;^\x,y). (1.10) 



In the following lemma we provide an upper estimate for the "massive" Green's 
function. It is well-known and easy to see that this Green's function decays exponen- 
tially in the distance to the origin. We prefer however to insert a proof for the sake 
of self-containedness, and indicate the power of 7 entering in the exponent. 

Lemma 1. There exist G > and c > such that for all A C Z*^, 7 < 1, and 
x,y eZ"^, xj^ y, 



G^;i\x,y)<G^^'\x,y) 
G^^\x,y)< 



y\ 


if \x 


- y\ 


> 


7 




if \x 


- y\ 


< 


7 




if \x 


- y\ 


< 


7 



c_ 

Ix-yl" 



Furthermore, there exists G' > 0, c' > such that the reverse inequalities hold with 
G' replacing G and c' replacing c. 

Proof. First note that G^j^\x, y) < G^^]{x, y), since in the Green's function G^j^\x, y), 
the random walk is killed upon exiting A. Next, we have 

G?i(x,y)= Y: [^) Pn{x,y)^ (1.12) 

n=\x—y\ 

where Pn{x,y) denotes the n-step transition probability of simple (nearest neighbor) 
random walk {Sn} on Z'^, and where the sum over n starts at n = |x — ?/| since the 
nearest neighbor random walk has to make at least that number of steps to reach y 
from X. 

We use that pn{x,y) satisfies the Gaussian upper and lower bounds: 

Pn(x,2/)<^e-^^l--^l'/", (1.13) 
and for |x — y| < n and n of the same parity as \x — y\, 

Pn(x,2/)>^e-^^|^-^|'/^ (1.14) 
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Since we do not know of a reference where these are stated in this form, we sketch 
the simple proofs. The upper bound follows from the large deviation bound PdS*™] > 
|2^|/2) < C3exp(— Ci|a;p/m), and the fact that ¥{Sm = y) l£ CimT'^^'^, a consequence 
of the local limit theorem. Taking m = [?t-/2J, the two imply: 

n—m n—m ) 

z:\z\>\x\/2 z:\z-x\>\x\/2 

< C^n-'^l'^ exp(-Ci|x|Vn). 

The Gaussian lower bound follows from a chaining argument: assuming |xp > n, let 
m = [\x\'^/n\ , and let i/q = 0,yi, . . . ,ym = x he points such that \yi-i — yi\ < 2n/\x\. 
Consider the balls Bi = B[yi, n/\x\). Then for any z G -Bj_i, the central limit theorem 
implies F{Sn/m ^ Bi\So = z) > c > 0. We get, using the local limit theorem, 

Pn{0,x) > d^-^CMV^y/^ > a.n-''/^ exp{-C[\x\yn) 

Inserting the upper estimate into fll.l2p . and using the notation C3 = I/7, we 
have 2d/ {2d + 7) < exp(— C37), we obtain 

71=1 ^ ' 

°° I 1 2 

^ C2j:^2^M-Cs^n-C,^-^). (1.15) 

n=l 

In the case |x| > 7"^''^, the bounds follow from estimating separately the sums 

Yl ;^^^P(-^il^lV^) and ;^^^P(-^37^)- 

l<n<\x\/^ n>\x\/^ 

(1.16) 

In the case |x| < 7 d = 3, the bound follows by ignoring the killing, and using 
[HI Theorem 1.5.4]. In the case d = 2 we again estimate the sums (11.161) . 

The proof of the lower bound is similar, starting from the lower bound on p„(x, y). 

□ 



1.3 Rational 7 

When 7 is rational, the dynamics can be reduced to that of a discrete model. Indeed, 
let 7 = k/n, with k > 0, n > 1 integers, and k and n relatively prime. Then the 
vector {nrix)x(^K changes by integer amounts both during addition of unit height, and 
during toppling. Hence the fractional part {{nrix})x<^A remains invariant, and can be 
"factored out". We define the map y^A : '^a — * {0, 1,2, .. . }^, by {(pi^{vi))x = [nrix]. 
The discrete toppling matrix is 

{2dn + k if X = y; 
-n if X ~ y; 

otherwise. 



6 



and it has associated addition operators a^^x-, allowed configurations 1Z\^ and top- 
pling operators T^^- It is easy to verify that {W,rjw) is a (/i;/n)-FSC if and only 
if {W,(p!^{riw)) is an FSC with respect to A. This implies that ip^iTl^^^"^) = T^a, 
and that the stationary measure for the discrete model coincides with (p^m^^^^'^ . The 
relation between the toppling operators is Ta,xV^a = ^hT^lf^ ■ Consequently, since 
adding unit height in the continuous model corresponds to adding n particles in the 
discrete model, with the notation 6A,a: = (^a.x)", we have ^a.^V^a = '/^aOax"^- 

The elements {6A,x}xeA generate the sandpile group for Aa- To see this, note that 
the order of the group, det(AA), is relatively prime to n. Hence, powers of 6a,x yield 
all powers of aA,a:; and the claim follows. 

1.4 Discretized heights, spanning trees 

The measure can be described in terms of discrete height variables. Consider 
the space 

fif " = {0, 1, . . . , 2rf - 1, 2d]^ 
We introduce the discretizing map 



V^A : ^^i"^ ^T"' 



'A 



defined by 



m if m < r]y < m + I, m = 0,1, ... ,2d — 1; 
2d ii2d<7]y<-f + 2d. 



We define ip : n^"''^ Vl'^^^^'^ analogously. We define TZj^^^^ to be the set of configura- 
tions 7] E ri^^'^'^ which are allowed. By Remark ([T]), for any i] E we have 

7/ E n^^^ if and only if ^^{ri) E TZf^'. (1.17) 

By a (7, A)-ce//, we will mean a subset of TZ^^^ of the form TZ^^^ fl for some 

^ E TZj^^'^^. It follows from the above discussion that is uniform on each cell, hence 
m^'' can be uniquely specified in terms of the measures of cells. Let u^"* := ifj^m^j^^ 

We proceed to give a description of ^ e TZf^'. Let 

N{0 = \{yeA:^y = 2d}\. 

Then 

z/i^^(0= (1-18) 
^ det(Ai'^)) ^ ^ 

which follows from the fact that under ipA^^^ discrete heights ^ {0, . . . , 2c? — 1} 
go to intervals of unit length, and heights C,x = 2(i to intervals of length 7, hence 
Volume (V'a^(O) = 7^^^^- 
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In order to study the infinite volume limit, we interpret (^) in terms of weighted 
spanning trees. For this we adapt to our setting the Majumdar-Dhar correspondence 
between allowed configurations and spanning trees [H]. 

Burning algorithm (Majumdar-Dhar construction). 
^a(?7). Set 

Wi = {y e K: iy <2d} = Wo\{y e K 
For t = 1, 2, . . . , we recursively define 

Wt+i = Wt\\yeWt:iy> 

The set removed from Wt to obtain Wt+i is called the set of sites "burnt" at time 
t + 1. In particular at time 1 the sites in A with height > 2d are burnt. 

By induction on t and (11.71) . no site in A\Wt can be contained in any FSC Hence 
we have n^oWt = if and only if ^ G TZf^' (if and only ii rj e 11^2^). 

Now consider the graph Z'' obtained by adding a new vertex uj to and connect- 
ing it to every vertex. Let us call the newly added edges special, and the rest of the 
edges ordinary. Now we identify all the vertices in Z*^ \ A with uo (and remove loops). 
Call this new graph A. In A, every ?/ G A is connected to uo by exactly one special 
edge. Additionally, boundary sites of A are connected to uo by one or more ordinary 
edges. We denote by -E(A) the set of edges of A. 

We define a spanning tree Ta(^) of A. First, for each y G Wq \ Wi (that is, when 
^y = 2d), include the special edge of y in the tree. Now for each y G Wi, let t{y) > 2 
be the index for which y G Wt-i \ Wt, and let t(a;) = 1. For y G Wt-i \ Wt, let 

r{y) = \{{z, y} G E{A) : {z, y} ordinary and t{z) = t — 1}|, 
n{y) = \{{z,y} G E{A) : {z,y} ordinary and t{z) <t}\. 

From the construction we have 

r{y) = r, n{y) = n if and only if 2d — n < rjy < 2d — n + r 

if and only if 2d — n < C,y < 2d — n + r. 

A one-to-one correspondence can be set up between the 2d directions of the ordinary 
edges and the values {0,1, . . . ,2d ~ 1}. This induces a one-to-one correspondence 
between the r{y) ordinary edges and the values {2d — n, . . . ,2d — n + r — 1}. Include 
in Ta(^) the ordinary edge {z,y} corresponding to the value of ^y. Since each vertex 
in Wt is connected to a unique vertex in Wt-i, T\{^) is a spanning tree. It follows by 
construction that the mapping 

Ta : e ^ Ta(0 



Fix ri G n)^', and let ^ = 



: ^y > 2d}. 
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is one-to-one and onto the set of spanning trees of A. 

Let /ij^^ denote the distribution on {0, under which a spanning tree t has 

weight 7^(*) / det(A^''), where N{t) is the number of special edges in t. By construc- 
tion, for each cell N{Ta{0) = and therefore, by (|rT8ll . 

det(Ai^O 

2 The stationary measure of the dissipative model 

We can view A as a weighted network, where ordinary edges have weight 1, and 
special edges have weight 7. Then /i^^ is the weighted spanning tree measure on 
this network, that is, the measure where the probability of a tree is proportional to 
the product of the weights of the edges it contains. Wilson's algorithm [TT] can be 
used to sample from this distribution. This is described as follows. For a path a, we 
denote by LE{a) its loop-erasure, that is the path obtained by removing loops from 
a chronologically. Consider the network random walk on A, that is, the reversible 
Markov chain that makes jumps with probabilities proportional to the weights. Let 
JFq = {u;}, and let Xi, . . . ,xx be an enumeration of the vertices in A. If J-'j-i has 
been defined, start a network random walk from Xj, and run it until it first hits J^j-i- 
Let (Tj be the path obtained, and let J-'j = J-'j-i U LE{aj). By Wilson's theorem, J^k 
is a tree with the stated distribution. 

Let 7 > 0. By the ideas in ||4|, Section 5], for any exhaustion Ai C A2 C . . . of Z'^, 
the weak limit 

lim /ii-^ =: (2.1) 

exists. We will be interested in sampling from /i'^'^^ when 7 > 0. For this, Wilson's 
algorithm can be used. When 7 = 0, > 3, this is the statement of |H Theorem 5.1]. 
The algorithm in this case uses simple random walk, and the "root" is at infinity. 
When 7 > 0, the algorithm uses the network random walk on Z*^, stopped when it 
hits uj. The tree generated by the algorithm then gives a sample from To see 
the convergence fl2.ll) . couple the algorithms in A„ and 7/, starting walks at a fixed 
finite number of sites Xi, . . . , x^. The network random walks can be coupled until the 
first time the boundary of A„ is hit. The coupling shows that the joint distribution 
of {LE{ai), . . . , LE{(Tk)) in A„ converges as — > cxd to the joint distribution in Z*^. 
Since this information determines all finite dimensional distributions of the spanning 
tree, the claim (12. ip is proven. 

Under /i*-'^'', there is a unique path vr^. in the tree from x to u. Similarly, there is 
a unique path ir^^x under Let us write Vx for the set of self-avoiding paths from 
X to in Z*^, and = {y '■ \y — x\ < 1}. /^From the correspondence Ta we obtain: 

Lemma 2. Under the map T^-^ , the height is a function of {'n'A y}y^j^^ only. 
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Lemma 3. For any 7 > 0, and any exhaustion Ai C A2 C . . . of TL^ , we have the 
unique weak limit 

lim z/^^ =: z/(^\ 



A„ 



Moreover, for all x G Z , there is a map : Hi/gA/'j: '^v ~^ {Oj ■ ■ ■ > such that 
the law of {^x}xei.'' under u^^^ is the same as the law of {hx{7Cy : y G J^x)}xeZ'i- 

Proof. Coupling Wilson's algorithm in A„ and Z'^, we get that the joint distribution 
of {'n'A„,xi, • • • , '^An,xk) converges to the joint distribution of {hxi, • • • , T^xk)- This and 
Lemma [2] shows that for any xi, . . . ,Xk, the joint distribution of {^xi, • • • , ^Xk) under 
i^/^^, converges to a limit, and it has the form stated. □ 

For the next lemma, we make m^'' into a measure on Q^'^'^ via the inclusion map 
i : n^^^ where 

T]^ if X G A; 
if x G Z'^ \ A. 



Lemma 4. For any 7 > 0, and any exhaustion Ai C A2 C . . . ofZ'^, the weak limit 

lim m^?^ =: m^"'^ 



"Ar, 

n— >oo 



exists, and we have m^'^\TZ^'^^) = 1. 

Proof. Recall the fact that m^'* is uniform on (7, A)-cells. Therefore, under the 
measure mj^'* and conditioned on the value of ^ = ipxiv) ^ ^^a'^'^'^' ^^e random variables 
{rjy — C,y}y(^A are (conditionally) independent, with distribution 

Tjy — iy ^ Uniform(0, 1) for y with Q < ^y <2d — 1; 
f]y ~ ^ Uniform(0, 7) for y with C,y = 2d. 

Now let C A. Then the above implies that if we condition on the values of {C,y}y&v 
only, then {rjy — ^y}yev still has the conditional joint distribution in (12.21) (under 
m^''). This, and Lemma imply that the joint law of {riy}y(zv under m^j converges 
weakly as n — > 00, in the space ^ly^ . In fact, the limit distribution also satisfies (12. 2p . 

This proves the weak convergence statement. Finally, the limit gives probability 
1 to the event {r]v ^ T^v^} for any finite V, because 

for all WdV. □ 



We now prove that for dissipation 7 > 0, we have exponential decay of correlations 
of local observables. In [TT] this was shown for large enough dissipation, in the discrete 
case. 

By a local function, we mean a function / : Q'^'"''' M (respectively / : Q^^^ R) 
such that / depends only on coordinates in a finite set A C Z*^. 
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Theorem 1. For all 7 > 0, there exist C = C^'^\c > 0, such that for all bounded 
local functions f, g : Q'^^^^^ — ^ M with dependence sets A, B, we have 

|E,w[/(7] -E,m[/]E,(,)[(7]| < C|A||5|||/|U||^?||ooexp(-Cv/^dist(A5)). (2.3) 

The same statement holds for the measure m^'^^ . 

Proof. First consider the case when / depends only on and g only on Then by 
Lemma El / and g depend on the paths {vr^j^g^^, and {TT^jioeATy, respectively. Use 
Wilson's algorithm starting with the vertices in and then using the vertices in 
f/y. We couple the random walks appearing in the algorithm to a new set of random 
walks {S'^}ui£j\fy as follows: 5*^ = Sw until the path hits {u} U (U^gA/^vr^), and it 
moves independently afterwards. Let {'rc'u,}wej\fy be the paths created by Wilson's 
algorithm started from Afy, using the S'^^s. Let g' be a copy of g that is the function 
of {tt^ : w G N'y}- Then the left hand side of (12. 3p equals ^[f{g — g')]. This is 
bounded by 2||/||oo||fi'||oo times the probability that g 7^ g'. The latter is bounded by 
the probability that the random walks used for g do not intersect any of the random 
walks used for /. This is bounded by 

E E (2.4) 

z&Afx w&Ny nSZ"^ 

where G'^^^'^^ is the Green's function for the discrete time walk that steps to each 
neighbor with probability 1/ (2(i + 7) and to uo with probability 7/(2^ + 7). Hence we 
obtain ([US]) from (fTTTD . 

In the general case, we can repeat the argument with the vertices 'm. M.f = U^jgAA/'x 
and }Ag = Uj^gsA/'j,. This leads to an estimate similar to (12. 4p . where now we sum 
over z E Aif and w G Aig. This implies the claim. 

Let us prove the statement for the continuous model. Due to the representation 
(D, if we set 

foiO = ^m(^)[f{vM{v)AuB = Uub] 

then we have E^(^) [f{v)g{v)] = Km [/o(0^o(0]> ^mh) [/(^)] = Kh) [/o(0] and E^(^) [g{ri)] = 
K(~i)[9o{0]- Now we can apply the discrete result to /o and go. □ 

Remark 2. TheoremUl shows that in the system with dissipation 7, the "correlation 
length" scales as l/^/y (as •y ^ 0). A natural length associated to the system with 
dissipation 7 is the size (diameter) of a typical wave. Using a two-component spanning 
tree representation, one can show that (see section 7 of JE/ for more details), starting 
from a configuration chosen from the stationary measure m^'^\ the probability that 
a typical wave contains x has the same asymptotic behavior (for large \x\) as the 
probability that a loop-erased random walk starting at x, and killed at rate 7 ever 
reaches 0. This probability equals G^'^\0, x)/G^'^\0, 0), which decays (up to polynomial 
factors) as e"'^'''^*-'^-' with ^(7) = Cj ^ by lemma{^ So we can think 0/1/^/7 as the 
"typical wave-size" as 7 ^ (or "avalanche- size" , which is the same, because the 
number of waves is of order 1). 
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The non-dissipative discrete sandpile measure v = u^^^ was constructed in [I] and 
[HI Appendix], as the weak hmit of ufK 

Proposition 2. We have hm^_»o z^*-"^-* = and hm^^o = fn^^\ 

Proof. For G ^ij^^'^'^ chosen from z/^^'*, consider the random field of maximal heights 
{hA,x}x&A ■= {I[^A,x = 2(i]}a;gA, and the random set Ha = {x e A : hA,x = !}■ Due 
to the correspondence Ta, f{x G i^A) equals the probability that the special edge 
containing x is included in Ta(0- Using Wilson's algorithm started at x and fll.lSp . 
we see that this probability vanishes in the limit 7 — > 0, uniformly in A. Hence for 
any finite V C Z'^, F{Ha (IV = 0) goes to 1 as 7 0, uniformly in A. 

By Remark [H after removal of the sites with C,a,x = 2(i, the joint distribution of 
the heights of the remaining sites is uniform on allowed configurations with heights in 
{0, 1, . . . , 2d — l}. Therefore, given Ha, the conditional distribution of {iA,x}x£A\HA is 
given by J^f^^^, the non-dissipative sandpile measure in A\Ha- Due to the convergence 

~^ ^^^^ ^ as W ^ Z"', for large V and on the event {Ha r\V = 0} the conditional 
distribution is close to z/'^°\ uniformly in Ha and A D V. The above observations 
imply that as 7 and A Z"', u^^^ — > z/*^°\ This implies the first statement of the 
Proposition. 

For the second statement, we again use the representation (12.21) . Then the conver- 
gence of finite-dimensional distributions follows from the first part of the Proposition. 
Tightness holds trivially since all the measures under consideration have support con- 
tained in (the same) compact set. □ 

Remark 3. The process {hx}x&'i fo-ct a determinantal process fT^ . that is there 
exists a kernel K{x, y), such that for n > 1 and distinct xi, . . . , x„ G Z'^, 

F{hx, = l,...,hx„ = l) = det{K{xi, Xj))l^^,. 

This follows from the Transfer Current Theorem ^ Section 4] applied to the collection 
of special edges. 



3 The dynamics of the dissipative model 

Let (f : Z"^ ^ (0,00) be a bounded function, and let {N^)t>o '■= {Afx,t}xezd,t>o be 
a collection of independent Poisson processes, where {N'^^^)t>o has rate v^(x). We 
assume these are defined on a probability space (X, JF, P) . We want to define the 
dynamics {rit}t>o of the dissipative model for an initial configuration 779 as the result 
of stabilizing tjq + Nf. That this is well-defined is the result of the next two lemmas. 

First we consider stabilization of infinite configurations. Fix 7] & X, and a sequence 
xi,X2, . . . of 7- legal topplings. 

Definition 1. A sequence xi, X2, ■ ■ ■ of'y-legal topplings of a configuration f] is called 
exhaustive, if for every n G N, and for every 'j -unstable site x of Tx„ o . . . o Tx^{ri), 
there exists m > n such that Xm = x. A sequence Xi,X2, ■ ■ ■ of j- legal topplings of a 
configuration t] is called stabilizing it the limit \imn^ooTx„ o . . . o Tx^{ri) is ''^-stable. 
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Note that every stabilizing sequence is exhaustive and if the configuration is 7- 
stabihzable, every exhaustive sequence is stabihzing. 

Just as in the finite case, the number of times each site topples in an exhaustive 
sequence (which may be infinite) is independent of the sequence (see [TJ Lemma 4.1] 
for a proof). It can also be seen the same way that if yi,y2, ■ ■ ■ is another 7-legal 
sequence of topplings, then each site topples at most as many times as in an exhaustive 
sequence. Call 

{N^'^^ri)x = the number of times x topples in an exhaustive sequence. 

We say that r] G X is 'y-stabilizable, if {N^^^rj)^ < 00 for all x G Z"'. In this case, we 
denote the stabilization by 

S^^\r]) = r]- A^^^N^^'>r]. 

Recall that {N^^ri)^ denotes the number of times x G A topples during stabilization 
in A. 

Lemma 5. Let 7 > 0, G A". 

(i)The function Nj^\ri) is monotone increasing in A and rj. 
(a) If rj is "y-stabilizable, we have N'''^^{ri)x = supj^ N^\ri)x, x G Z'^. 
(Hi) If 7] is 'f-stabilizable, we have S^^\ri) = lim^ S^^\ri) G fi*^'^^. 
(iv) If (1,(2 € , o-nd C^i + C2 is 'y-stabilizable, then we have 

5(^HCi + C2) = 5W('5("HCi) + C2). (3.1) 

Proof, (i) Consider A C A'. Stabilize first rj in A, and record the amount of height 
fiowing out of A. Now stabilize in A', with further topplings, if necessary. To prove 
the other statement, consider rj < rj'. Stabilize first rj. Adding height rj' — rj does not 
affect the legality of the sequence of topplings. Hence we can stabilize rj' by further 
topplings if necessary. 

(ii) Fix a sequence Ai C A2 C ... such that UfcAfc = Z'^. Stabilize 77 in Ai, 
then in A2, and so on. We thus get an exhaustive sequence, and each site x topples 
sup^(A^|'^^?7)^ times. 

(iii) Note that if r] is 7-stabilizable, we have N'^'^'^rj = limx{N^/ij'' t]) , and hence by 

S^^\r]) = r]- A^^^N'^^'^ri = lim 

(iv) This is essentially not more than the Abelian property. However, since there 
can be infinitely many topplings in stabilizing (i, care needs to be taken. Note that 
Ci < Ci + C2 implies that any legal sequence for (i is legal for (i + (2 as well. In 
particular, ^1 is also 7-stabilizable. Fix a stabilizing sequence 1/1,1/2, • • • for d. We 
now construct an exhaustive sequence for Ci + C2- Fix again Ai C A2 C . . . . We write 

d+A = {x G A'^ : 3y G A such that x ~ y}; 
d-A = {y e A : 3x e A"" such that x ~ y}; 
A = A U d^A. 



lim5i^)(r/). 
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Let mi be an index such that for all ?/ G Ai 

|{l<A;<mi:y, = |/}| = (iVWCi)r 

Note that . . . , Vmi is 7-legal for Ci + (2- Now stabilize Ty^^ o . . . Ty^{Ci + C2) in Ai \ 
d-Ai, recording the amount of height flowing out of Ai \(9_Ai. This leads to potential 
further topplings in Ai \ d^Ai at 2:1, . . . , z^^. By construction, yi, . . . , Um^, zi, . . . ,Zn-^ 
is 7-legal for ^1 + and Zi, . . . , z^^ is 7-legal for S^^^{C,\) -|- C2- Moreover, the con- 
figuration in Ai \ (9_Ai after topplings at . . . , j/mi, Z\, . . . ,Zn^ coincides with the 
stabilization of >S(^)(Ci) + C2 in Ai \ 9_Ai. 

Now select an index > vtlx such that for all y € A2 

|{l<A;<m2:yfc = y}| = (ArWCi)r 
Since the topplings at 2:1, . . . , do not change the configuration in A^, the sequence 

Dl-i • • • -I Vmi i ^li ■ ■ ■ 1 ^ni ; Vmi + l ; • • • ; l/m2 

is 7-legal for Ci + C2- Now we stabilize in A2 \ 9_A2 via topplings at Zm+i, . . . ,Zn2- 
Again by construction, 

is 7-legal for Ci + C25 and 

■^1 ; • • • ; ) + 1 ; • • • ; ^n2 

is 7-legal for S^^'^Ci) + (2- Continuing this argument, we obtain an interlacement of 
two sequences yi, y2, ■ ■ ■ and zi, Z2, . . . which is exhaustive for ^1 + ^2, and hence its 
final result is S^"'\Ci + (2)- On the other hand, the configuration in A^ \ d^A^ after 
the topplings at 

yii ■ ■ ■ 1 Vmii Zl, . . . , Znn ■ ■ ■ ?/mfc_i+l) • • • ; Vmi^i -^n^.i+l; • • • ; ^rik 

coincides with the stabilization of ~^ C2 in Afc\9_Afc. This yields the claim. □ 

Lemma 6. Let 7 > 0. 

(i) If 1] E X satisfies Yliy&'i G^'^\x,y)r]y < 00 for all x E Tf- , then i] is j-stabilizable. 
(a) With F-probability 1, for any 7]q E fi'-'^-' and any t >0, r]Q + Nf is 'j-stabilizable. 



Proof (i) By (O, 

A^;^^Nj^^ = r]-S'j^\r]) <r]. 

Hence, by f fTTT]) 

sup(A^|'^^?7)^ < sup ^ g57V> < '^^^\^^y)^y < (3.2) 
^ ^ 2/eA jyez'* 

(ii) Due to the boundedness of (p and estimate fll.lip . the configuration rjQ + A^f 
satisfies the condition in part (i) for all t > 0, with probability 1. □ 



14 



Due to Lemmas and El the process 

Vt = S^''\vo + Nn. (3.3) 

is well-defined for any initial configuration rjo E The computation in (13.21) also 
gives that the addition operators 



a^r/ := S^^^t] + 5,) = lim4'^^(r/ + 5,.) 



are defined for any rj G 

Recall that n^-* (x, y, rjo) denotes the number of topplings occurring at y in comput- 
ing S^^^ (?7o + We also define n^'^'' (x, y, rjo) := sup^ n^-* (x, y, rjo) and the avalanche 
started at x (in A and in Z*^) by 



AvaU(^) = {?/ e a : n^^(x,i/,r7) > 1} 
Av(?')(r/) = {y G Z'^ : ^(^^(x, r]) > 1}. 

Lemma 7. Let 7 > 0. 

(i) We have m(^) (lAv^'^^l < 00) = 1. 

(a) The transformation ai^^ leaves m^'^'^ invariant. 

Proof. Take r/o distributed according to m'^'^\ Due to (11.91) . 

(x, ?7o)] = lip E^(^) [n J'^ (x, y, r^o)] 

A 

= limlimE M [n!^'^^(x, 

A V "V 

- ^i^^m(7'[^?^(^'2/'^o)] 
= G^^^\x,y). 



(3.4) 



(3.5) 



Inequality (13. 5p and the estimates (II. lip . (ll.lOp imply E^(7) |Av^'''-'(?7)| < 00, for any 
X, implying (i). By a simple argument [SI Section 4], (i) implies (ii). □ 

Theorem 2. Let 7 > 0. 

(i) The process {rjt\t>Q is Markovian. 

(ii) With F-probability 1, for any r]o G ^l^^^ we have 

r]t = lim JJ (ai"^^)^"'* (vo)- 

x£V 

(Hi) If rjQ is distributed according to m^"''^ then {'r]t}t>o is stationary. 

Proof, (i) By (13. 3p . for < t < t + s, and rjo G rjt and rjt+s are well-defined, and 
Tjo + Nf, r]Q + Nf^g are 7-stabilizable by Lemma [U] (ii). By Lemma [S](iv), we have 

r^i+. = 5W(r/o + Nf^;) = S^^\vt + N^,, - Nf), s > 0. 
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This implies the Markov property. 

(ii) Condition on a reahzation of the Poisson processes such that for all y & 
'^x&zd ^x,tG^^K^i y) < c>o. Let xi, X2, ... be an enumeration of Z'^, and let 

n 

Vn = VO + ^N^^/xi, 
i=l 

r] = r]o + Nf] 

n 
i=l 

Note that C is indeed well-defined by Lemma EJ^ii). Given W O Z'^ finite, select A 
such that {Nj^^ri)y = {N'''^'^ri)y for all y G W. If n > no (A), we have //„ = in A, and 
therefore 

where the first inequality is due to Lemma [5]^ii), and the second to Lemma [5t^i). This 
shows that {Cn)y = Cy fo^ U ^ ^) wheu n > rio(A). Hence C,n proving (ii). 

(iii) Due to Lemma [Tl^ii), and part (ii), rjt is an almost sure limit of configurations 
distributed according to m*-^-*. This completes the proof. □ 

We now consider the dynamics in the case 7 = 0. For the discrete model, this was 
constructed in [8j. Due to the observations in Section [L3l the discrete construction 
immediately implies the following statements for the continuous model. 

Theorem 3. Let 7 = and d > 3. 

(i) For m^^^-a.e. rjo, the configuration % + 5^ is 0-stabilizable for all x E TL'^ . 

(ii) a>x^ leaves mP^^ invariant for all x G Z''. 

(iii) Assume that satisfies 'YlixeZ'^ (p{x)G^'^\x, 0) < 00. Then m^^^ ®P-a.s., the limit 

T]t = lim JJ (4°))^"''r7o 

x£V 

exists, and equals S^'^\r]Q + N^'). 

We denote m := m^^\ G(-, ■) := G'(°)(-, ■) and Av^iv) ■= Avf (r/). 

4 Zero dissipation limit of the stationary processes 

Let (fi be an addition rate such that 

J2v{x)G{x,y) <oo (4.1) 
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for all ?/ G Z"'. Let rjl denote the stationary process obtained from starting from 
Vq''' •= v'^'^^ distributed according to m'^'^\ addition according to independent Poisson 
processes with rate ^{x) at x G Z*^, and stabilizing with dissipation 7. Similarly, 
let rjt denote the process starting from rjo = r] distributed according to m, addition 
according to independent Poisson processes with rate (p{x) at x G Z'^, and stabilizing 
without dissipation, i.e., with 7 = 0. 

Theorem 4. The process ri[^^ converges weakly in path space to rjt. 



We need the following two lemmas. 
Lemma 8. Suppose that < 7' < 7. Then 

(zi) Avo-f{r]) C Avo(?7). 



Proof, (i) Consider an exhaustive sequence yi,y2, ■ ■ ■ of 7-legal topplings for rj. We 
show that the same sequence is 7'-legal for rj. Since yi is 7-unstable in rj, it is also 
7'-unstable. After its 7-toppling let us add height 7 — 7' > at This has the 
same effect as if we performed a 7'-toppling, and the added extra height does not 
affect 7-legality of the sequence. Adding similarly after each 7-toppling shows that 
yi,y2, • • • is 7'-legal, and (i) follows by the remarks preceding Lemma El 

(ii) This follows immediately from (i). □ 
Lemma 9. Suppose that 7„ | 0. Then for m-a.e. rj, 

Mm a'^^'"\7]) = a^o\v)- (4-2) 



n— »oo 



Proof. We have 

m(|Avi°)(77)|<oo) = l; (4.3) 

see [H Theorem 3.11]. Also, m{rix ^ {0,1, . . . ,2d — 1}, x E Z'^) = 1. On the intersec- 
tion of the two events we have N^'^"\r]) N^^^(r]), which together with A*^'^"^ A^^^ 
implies the statement. □ 

Proof of Theorem^ Recall the definition fll.Sp of the metric on [0, 2d -\- 7)^'', that is 

dist(r7,C)= ^2^-lmin{|r/,-C.|,l} 

As a first step, for a given 5 > 0, we will prove that for all e > 0, there exists a 
coupling Ai^'^'^ of nS'^'^ and m such that, in this coupling, with probability at least 

dist(4'^)(r^i),4^)(r^2))<5 (4.4) 

This will allow us to deal with the convergence of the processes for addition rates 
with finite support. 
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First, for given 5, e > 0, we choose V large enough such that if rj, ( agree on V, 
then dist(r7, < 6, and such that 

m (|A^| ^V'^<e (4.5) 

Such a choice of V is possible, since avalanches are finite with m-probability one 
by fl4.3p . By Lemma [H](ii), we then have the same estimate (14.51) for Avq^''(?7). 

Since m^'^^ — > m weakly, and restrictions to finite volumes A C Z*^ of m^'^^ and m 
are absolutely continuous with respect to Lebesgue measure, by [15], Proposition 1] 
we have the existence of 70 > such that for all 7 < 70 there exists a coupling A^'^'^^ 
of m^'*'^ and m such that 

M^^\vl = vl, e V) > 1 - e (4.6) 
In the coupling Ai^"'^ we then have, by fl4.5p 

for aW y G V, with probability at least 1 — 3e. Therefore, the probability that the 
distance dist{ap\'ri'^) , a^Q\'ri'^)) is larger that 6 is smaller than 3e. 

So far, we can conclude that 

4^)(r^W)-.4^)(r^°) 
weakly as 7 — 0. By Lemma [HI al^\ri) — » aQ{ri) for m-a.e. r]. Hence we have 

a^^\v^''^) ^ aoiv') (4.7) 

weakly as 7 — >• 0. 

Analogously, using finiteness of avalanches, we conclude that for any finite set 
B C Z"^, and natural numbers Ux, x G -B, we have 

n(4^))"^(r/(^))-n<^(^°) (4-8) 

weakly, as 7 ^ 0. 

Therefore, we have convergence of the processes rjl"'^ rjt for addition rates with 
finite support, i.e., such that '^{x) = for x ^ D with D C 'Z'^ finite. 

The next step is to pass to general addition rates; we use the convergence of 
semigroups argument, as in the proof of [11], Proposition 4.1]. 

Let Sf""^ denote the semigroup of the process 77^^'' with addition rate ip, and Sf 
the semigroup of the process rjt (with zero dissipation) with addition rate ip. Both 
semigroups are well-defined as long as p has finite support. 
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For a local function /, with dependence set Dj, and for addition rates ip,ip' of 
finite support, we have the estimate (51) in the proof of [llj Proposition 4.1] 

- Sf'\f)\ < Ct J2J2G^^\x,y)Hy) - ^'iy)\ 

xeDj y 

x&TTf y 

and 

^m\St{f) - Sf{f)\ < Ct ^5^G(x,y)|v.(i/) - v\y)\ (4.9) 

Note that this estimate in [H] is given in the context of a model with discrete heights. 
However, it is based only on the estimate f ll.lOp for the numbers of topplings, which 
is valid for the continuous height model, and therefore it extends directly to the 
continuous height model. 

Hence, if for a sequence (y?'-"^ of addition rates of finite support, for an addition 
rate if (not necessarily of finite support) and for all x ^if" 



5^G(x,i/)|y.(y)-^(")(l/)H0 



as n ^ oo, we conclude for all local /, 7 > 0, Sf '^{f) is a Cauchy sequence in 
L^(m^^^), and hence converges to \I' := Sf'^{f) for all 7 > 0, as n — ^ 00. This 
semigroup Sf'^ then defines a corresponding stationary Markov process r]f'^'^\ and rjf 
(for 7 = 0). 

As the convergence of the semigroups implies the weak convergence of the corre- 
sponding stationary processes, we conclude for all 7 > 0, 



and 



vf'' - Vr (4.10) 



vf'^vf (4.11) 



weakly on path space as n — 00. 

Therefore, if satisfies (14. ip . let ip^"^^ denote Lp^'^\x) = ip{x)I{x E [—n,nY), then 
we have for all n G N 

Vf Vi (4.12) 

as 7 — > 0. Combination of (14 . 1 p . (14 . 1 1 1) . (14 . 1 21) together with a three epsilon argument 
then concludes the proof. □ 
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